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We model atomic motion in a sliding superlattice potential to explore topological “charge pump¬ 
ing” and to find optimal parameters for experimental observation of this phenomenon. We analyti¬ 
cally study the band-structure, finding how the Wannier states evolve as two sinusoidal lattices are 
moved relative to one-another, and relate this evolution to the center of mass motion of an atomic 
cloud. We pay particular attention to counterintuitive or anomalous regimes, such as when the 
atomic motion is opposite to that of the lattice. 
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Introduction - Slow periodic changes in a lattice po¬ 
tential can transport charge. For a filled band, the in¬ 
tegrated particle current per cycle in such an adiabatic 
pump is quantized [1]. We study a simple but rich ex¬ 
ample of this phenomenon, namely charge transport in 
a sliding superlattice, and draw attention to its coun¬ 
terintuitive properties such as regimes where the charge 
moves faster than the potential, or even travels in the op¬ 
posite direction. We argue that this anomalous transport 
is observable in a cold atom experiment. 

The quantum mechanics of particles in a one¬ 
dimensional (ID) superlattice is rich, displaying a frac¬ 
tal energy spectrum [2], and for incommensurate peri¬ 
ods boasting a localization transition similar to what 
is seen in disordered lattices [3]. While recent studies 
have focused on the tight-binding limit (the Aubry-Andre 
model) [4-21], we study the continuous limit of the ID 
superlattice, where, because of the weak potential, the 
single-particle spectra can be calculated perturbatively. 
Related cold atom proposals on quantized transport [22- 
26] have focused on the simplest superlattice where one 
sub-lattice constant is half of the other, and are therefore 
not in the anomalous regime which interests us. 

The ID superlattice can be mapped onto the Harper- 
Hofstadter model [2, 27]. The topological numbers 
(Chern numbers) associated with charge pumping can be 
mapped onto quantized Hall conductances [28, 29]. Re¬ 
cent experiments involving artificial gauge fields on 2D 
optical lattice have aimed to measure these 2D Chern 
numbers [30-34] . There are also related studies based on 
measurement of Hall drift [35], Bloch oscillations [36, 37], 
Zak phase [38-40], time-of-flight images [41-43], edge 
states [44-49], or density plateaus [50, 51]. 

In this letter, we study the charge transport in a ID 
sliding superlattice, where the moving lattice period is 
an arbitrary rational multiple of the static lattice. We 
analytically calculate energy band gaps and the topolog¬ 
ical invariants which give the integrated adiabatic current 
per pumping cycle [1]. The fact that this current can be 
made arbitrarily large and/or opposite to the direction 
of the sliding potential is counterintuitive. We present 
a physical interpretation of this phenomenon in terms 


of the quantum tunneling of Wannier functions between 
minima in the potential. We propose an experiment to 
detect this anomalous adiabatic current, and derive the 
optimal parameters. Through numerical simulations, we 
confirm that a negative integrated current and a non¬ 
trivial Chern number C = — 1 is readily measured in an 
experiment. 

Model - We consider the Hamiltonian of a ID super¬ 
lattice where one lattice adiabatically slides relative to 
the other, 

H = jdx^\x) + Vi{x,ip) + V 2 {xfj i){x) ( 1 ) 

where 2 p{x) represents the field operator of the particle, 
h is Planck’s constant, and m is the mass of the particle. 
The periodic potentials Vi{x,(p) = 2vico^{px — (p) and 
V 2 {x) = 2 v 2 Cos{qx) are commensurate, with lattice con¬ 
stants 271/p and 27r/q, intensities vi and V 2 . We take the 
relative phase p to be slowly varying in time. The period 
of the Hamiltonian is set by the greatest common divisor 
of p and g, i.e., n = gcd(p, g), as illustrated in the inset 
of Fig. 1. Treating 1/k as the unit length, we redefine 
XK ^ X, p/ K ^ p, and q/ k ^ q. Treating /2m 

as the unit energy, we redefine H/Er H, v\/E^ 
and V 2 /Er V 2 . The dimensionless Hamiltonian in the 
momentum space is then 

j^2 

H = Y y V’F* + (i'ie“*‘^V’Ffc+p + + h.c}j 

k 

( 2 ) 

Here fjk = f dx with dimensionless system 

length L and dimensionless momentum k. Since states of 
momentum k are coupled only to those of momentum k + 
n for integer n, we restrict ourselves to the first Brillouin 
zone {0 < k < 1) and rewrite the Hamiltonian, 

oo ^ 

0<k<l n=—oo 

+ {vie~^'^'tpl'tpn+p + V2pn'>Pn+q + h.C.) (3) 
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FIG. 1: Band structure of a ID superlattice for p = 2, g = 3, 
showing energy E vs dimensionless wave-vector k for weak 
potentials. Inset shows the two potential making up the su¬ 
perlattice, and illustrates the unit cell with period set by the 
greatest common divisor n = gcd{p,q). For this choice of p 
and g, the energy gap between the third and fourth band is 
set by the potential strength 2v2, the gap between the sec¬ 
ond and third band is set by the potential strength 2vi , and 
the small gap between the second and third band scales as 

~ ViV2. 


where we have suppressed the k index, writing jJjn = 
V^/c-l-n- 

To illustrate the resulting band structure, we impose 
a cut-off on n, and numerically diagonalize the Hamilto¬ 
nian in Eq. (3) for p = 2 and q = 3. The lowest four 
energy bands are shown in Fig. 1, and even for this sim¬ 
ple case the gaps display a range of behaviors for small 
vi and V 2 . The gap between the third and fourth band 
is induced by the potential V 2 (x), and is proportional to 
V 2 for weak potentials. The gap between the second and 
third band is induced by Vi(x), and is proportional to 
vi. The small gap between the second and third band 
is induced by the combination of these two potentials, 
which scales as ^ viV 2 . In the following section, we will 
discuss the origin of these scalings in the context of un¬ 
derstanding the lowest energy gap. 

Band gaps and topology - The eigenstates of the Hamil¬ 
tonian in Eq. (3) can be found perturbatively in the 
limit of ^ 1,^2 ^ 1. Suppressing the index /c, we write 
H = Ho d- Ai^i, with 

oo ^ 

Ho = + ( 4 ) 

n= —oo 

XHi = XpHp + X-pH-p + XqHq + X-qH-q, ( 5 ) 

where Hp = Er=-oo V’iV’n+p, Hg = J2n=-oo^i'^n+q, 
and A is a formal small parameter, with Xp = X'^p = 
and Xq = Al^ = V 2 . 

Eor small A and 0 < /c < 1, the eigenstates of the 
lowest band will be superpositions of | — 1) and |0), where 
\m ) =i’ln I vac). We let 5k = k — 1/2 denote the distance 
of k from the band crossing point and assume Sk > 0. 
The physics for 6k < 0 is analogous. While ordinary 


perturbation theory works far from the crossing {6k ^ e, 
where e will be precisely defined below), one must use 
higher order degenerate perturbation theory to find the 
eigenstates for Sk ^ e. As argued in our supplemental 
material [52] , the resulting effective Hamiltonian is of the 
form 


H,s = PHoP + E ( 6 ) 

>0 

r_|_,r_ >0 


where P = \ — 1)(—1| + | 0 )( 0|5 and s+,s_,r+,r_ are in¬ 
tegers. The operator is the contribution 

to i^eff involving the absorption of 77 = sp rq units 
of momentum from the lattices, where s = S-^ — S- 
and r = — r_. By conservation of momentum, 

a = (—= 0 unless 77 = 1. We linearize 
i^eff about Sk = 0, and write the operators in the basis 
{| — 1), |0)}. At the lowest nontrivial order, we have 


^eff 


— ^Sk aAe^^ 
aAe~^^ ^Sk 


+ const.. 


( 7 ) 


where A = and 5 ^, Vm correspond 

to the absolutely smallest solution to the Diophantine 
equation sp -\- rq = 1. This result agrees with a similar 
perturbative analysis carried out by Thouless et. al. [28] 
for a related model. 

The off-diagonal terms of Eq. (7) split the energy de¬ 
generacy at Sk = 0 , and create an energy gap of size 
AEg = 2\aA\. Eor example, if p = 2 ,g = 3, the abso¬ 
lutely smallest solution to the Diophantine equation has 
Sm = = I 5 as —p q = 1. Thus the energy gap 

is 2\aviV2\, as denoted in Eig. 1. Eor larger |s^| and 
|r^|, the energy gap can be extremely small. Ordinary 
perturbation theory would have sufficed in the regime 
where Sk ^ 2 |(aA|, allowing us to identify e as 2\aA\. 
Properties of higher bands can be analyzed similarly. 

By analyzing Eq. (7), we find that the lowest energy 
eigenstate of Eq. (4-5) has the form 


= —sin—— 


\k,(p) = -sin-e' 


1 ) + cos-e 


P --ix /2 


| 0 ) + ..., ( 8 ) 


where tan/3 = —2aA/Sk. The neglected terms are higher 
order in vi and V 2 . Eor Sk ^ 2 |(aA|, sin^ 1 and 
cos^ ^ 0 , and the coefficients are featureless. 

Slowly changing ip generates an adiabatic current [ 1 ]. 
Eor a completely filled band, the integrated current in 
one pumping period {p from 0 to 27r) is [53] 

p27T 

AQ = 2ttC= dk dpQkcp (9) 

Jo Jo 

where the Berry curvature is 


= i{d^{k,p\dk\k,p) - h.c.) = ^dkCos/S. ( 10 ) 
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We see is concentrated near the location of the en¬ 
ergy gap. Integrating the Berry curvature is trivial, yield¬ 
ing the Chern number C = Sm- Although our argument 
requires that vi and V 2 are small, due to the quantized 
nature of C, the result should hold for all nonzero vi and 
V 2 . In our numerical calculations with larger vi,V 2 , we 
find the curvature is roughly uniform over the Brillouin 
zone, but as expected its integral is unchanged. 

Anomalous charge pumping - By appropriately choos¬ 
ing p and one can make C = 5 ^ an arbitrary in¬ 
teger [54-57]. This means that in one pumping cycle 
a single particle may move arbitrarily far and/or oppo¬ 
site to the direction of the sliding potential. Such long¬ 
distance and/or retrograde transport seems unphysical. 
The magic comes from the adiabatic process: If the po¬ 
tential moves sufficiently slowly, the particles always stay 
in a global minimum of the potential. Due to the struc¬ 
ture of the superlattice, a slight motion of the potential 
could result in a dramatic change of the locations of the 
global minima (see Fig. 2(a)). Within a small portion of 
a pumping cycle, the particles may “tunnel” to the new 
global minima which could be a large distance away from 
the old minima. 


Viix,p) + V2ix] 



1 x/2tt 2 


FIG. 2: (Color online) (a) Illustration of adiabatic charge 
transport in a ID superlattice, where the particle “travels” 
through three unit cells to the left when the lattice poten¬ 
tial moves to the right by one period. Solid lines show the 
potential Vi{x,(f) -h V 2 {x) for different values of p. Arrows 
schematically show how the locations of the minima shift dis- 
continuously. (b) illustrates evolution of two separated poten¬ 
tials of the superlattice: the right-sliding potential Vi(x,p) 
(solid red) and the static potential V 2 {x) (dashed black), (c) 
Evolution of Wannier function. Arrows indicate the “tunnel¬ 
ing” process, (d) shows the evolution of integrated adiabatic 
current as a function oi p. In these plots we choose p — 2 
and ^ = 7, so the Chern number is C = Sm = —3. Other 
parameters are vi = 0.5 and V 2 — 0.25. 


To further quantify our interpretations, we calculate 


the integrated current 


A(^((^) — ^ dk j dp ^ 


and the Wannier function at lattice site j 


0</c<l 


where the Bloch wave function is 




1 

7z 


00 

n= —00 


( 11 ) 


( 12 ) 


(13) 


Here we choose a smooth gauge for the Bloch wave func¬ 
tion, so the Wannier function is well localized [58]. 

Fig. 2(d) shows the integrated current as a function 
of p^ calculated from Eq. (11) using a similar method 
to Ref. [59]. We see the function is “step-like”: Flat 
regions correspond to slow transport, while the particle 
motion is rapid in the steep regions. This is further illus¬ 
trated by the Wannier function in Fig. 2(c). During the 
slow transport, the Wannier function slowly drifts, while 
during the rapid transport, one peak drops in amplitude, 
and a second peak rises. This corresponds to tunneling. 

For small vi^V 2 ^ the timescale for adiabaticity r is re¬ 
lated to the size of the gap, 1/r |q;A| 

Thus when the Chern number C = 5 ^ is large and 
the potentials are weak, adiabaticity is hard to main¬ 
tain in a practical experiment. For large vi^V 2 ^ the gap 
again falls, owing to the large potential barriers. Fig. 3 
shows the energy gap as a function of vi and V 2 

for p = 2,g = 3. The gap has a maximium value of 
/S.Eg ^ 0.09 at vi = 0.23 and V 2 = 0.95. An optimized 
experiment would be performed with these parameters. 



FIG. 3: Energy gap AEg as a function of I’l, i ’2 for p — 2^q — 
3. The gap has a maximium value of AEg 0.09 at r'l =0.23 
and V 2 — 0.95. 
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Experimental proposal - To observe this anomalous 
current, we envision a Fermi gas confined to a quasi- 
ID tube, such that only one transverse mode is occu¬ 
pied. Although the present analysis is ID, we expect 
the phenomena will persist for more general transverse 
confinement. Along the tube we engineer two longitu¬ 
dinal periodic potentials Vi(x, ip) = 2vi cos{px — ip) and 
V 2 {x) = 2 v 2 cos{qx) via two pairs of counter-propagating 
laser beams. The time-dependent phase ip = Suj t is pro¬ 
duced by a frequency difference Suj between two of the 
beams. To satisfy the adiabatic condition, we require 
hduj The resulting adiabatic particle current 

can be detected by observing the motion of the center of 
mass of the cloud: After time t = 2t:N/5uj^ the center 
of mass should move a distance Tc = 2'kCN j k. A dimen¬ 
sionless measure of this displacement is Xc — K,rc. The 
displacement can be measured in-situ [60-62] or after 
time-of-flight [63]. 

In modeling this experiment, one must account for the 
finite cloud size. We include this physics by adding a 
harmonic potential along the tube, V{x) = muJ^x^I2. 
Such potentials are always found in such experiments. 
Within a local density approximation, the lowest band 
will be filled at the center of the cloud, but only partially 
filled near the edge. Although our Chern number argu¬ 
ment only applies to the central region, we still expect 
the center of mass motion to be nearly quantized. For 
ficjo ^ and the particle number is much greater 

than one, only a very small portion of particles live at 
boundaries. Our numerical simulations (detailed below) 
confirm this resuits. For a typicai experiment, ujq ^ 10 
Hz, and vi/h^V 2 /h 100 kHz [64]. 

Because of the trap, the dispiacement Tc cannot be 
made arbitrarily large. When mu{^r\j2 is of order of the 
band gap atoms can tunnel to the higher bands. In 

our numerical simulation, we see that for small the 
maximum displacement scales as 1/ujq. 

Numerieal simulation - In order to see the feasibil¬ 
ity of our experimental proposal, we numerically sim¬ 
ulate the dynamical evolution of a ID Fermi gas. We 
take the many-body state to be a Slater determinant, 
made up from single-particle wave functions '0i(x, t) with 
1 < i < z/, where u is the number of fermions. At time 
t = 0, 'ipi{x,0) is the ith eigenstate of the Hamiltonian. 
We evolve %jji{x^t) via the time-dependent single-particle 
Schrodinger equation, and then calculate the center of 
mass Xc{t) = x\^lji{x^t)\'^dx. Fig. 4 shows 

the results for p = 2, g = 3 where the Chern number 
is C = —1. We see Xc < 0, meaning that the parti¬ 
cles travel in the opposite direction to the sliding poten¬ 
tial. Remarkably this retrograde motion persists even 
for relatively large 6uj. As ^ 0 the motion becomes 
quantized. A typical experiment has Er/h ^ 100 kHz 
[64], so the Chern number C = — 1 is readily extracted 
when Suj < 200 Hz. The inset of Fig. 4 shows the evo¬ 
lution of the center of mass in one pumping cycle for 


Xc/27r 



FIG. 4: (Color online) Displacement of the center of mass (in 
units of the super lattice) after one pumping period T — 2'kISuj 
for z/ = 63 fermions in a superlattice with p — 2^q — — 

0.23Fr,r’2 = 0.95Fr, and a harmonic trap/icjo = 2.2xlO“^E'r. 
Physically, Suj is the detuning between the beams producing 
the lattice with wave-number p. We see Xc/27r ^ C — —1 as 
Suj decreases. Inset shows the evolution of the center of mass 
for hSuj = 0.002F;^. [c.f. Fig. 2(d)] 

hSuj = 0.002F^^. We see the function is “step-like”, simi¬ 
lar to the ideal case (no harmonic trap and adiabatic) in 
Fig. 2(d). 
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SUPPLEMENTAL MATERIAL 

Here we derive an effective Hamiltonian for Eq. (4-5). 
For small A Since the eigenstates of the lowest band will 
be superpositions of | — 1) and |0), motivating projection 
operators 

P = |-1)(-1| + |0)(0| (14) 

Q = 1-P. (15) 

The states \m) = T/^^l^lvac), satisfy H()\m) = ^{kx ^ 
m)^|m). We seek eigenstates H\ 2 p) = E\ 2 p). We break 
the wave function into two parts 

\2p) = P\2p) + Q\i;) = \2po) + IV'ex), (16) 

where \ 2 po) is in the low energy sector, and \^ex) is a 
superposition of the higher-energy states. The eigen- 
equation is then decoupled into two equations 

PH\i;) = PE\^) = E\^o) (17) 

QH\^) = QE\^)=E\^,^). (18) 
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Inserting the identity = P + Q = 1 on the left 

hand side of Eq. (17)-(18) and substituting j-^ex) in terms 
of we obtain a closed equation for It/^o), 

^efflV'o) = (19) 

where 


Using the identity PHqQ = 0 and expanding the second 
term of Eq. (20), we obtain 


H,s = PHP + PHQ^-^j^QHP. (20) 


H.. =P/f„P + APH.P + A“PH.o|: (21) 

j=0 


This equation can be written as 

H,ft = PHoP+ 

S+,S_ >0 

r+,r_ >0 

( 22 ) 

where the momentum conservation implies that a = 
(—|0) = 0 unless sprq = 1, where 
s = — 5_ and r = — r_. In our problem, the lowest 

order contribution to a has either 5+ = 0 or s_ = 0. 
Similarly r+ = 0 or r_ =0. The lowest order contribu¬ 


tion to the diagonal elements of Peff corresponds to an 
identity matrix. 

Linearizing Peff about Sk = 0^ and writing the opera¬ 
tors in the basis {| — 1), |0)}, we have 

= ( aTeSx ) + const-’ (23) 

where A = , x = and correspond 

to the absolutely smallest solution to the Diophantine 
equation sp rq = 1. 







